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Abstract 

Extending results of Oh and Zumbrun in dimensions d > 3, we establish nonlin- 
ear stability and asymptotic behavior of spatially-periodic traveling-wave solutions of 
viscous systems of conservation laws in critical dimensions d = 1,2, under a natural 
set of spectral stability assumptions introduced by Schneider in the setting of reaction 
diffusion equations. The key new steps in the analysis beyond that in dimensions d > 3 
are a refined Green function estimate separating off translation as the slowest decaying 
linear mode and a novel scheme for detecting cancellation at the level of the nonlinear 
iteration in the Duhamel representation of a modulated periodic wave. 

1 Introduction 

Nonclassical viscous conservation laws arising in multiphase fluid and solid mechanics ex- 
hibit a rich variety of traveling wave phenomena, including homoclinic (pulse-type) and 
periodic solutions along with the standard heteroclinic (shock, or front-type) solutions 
\GZ\ \Z6\ \OZl\ IOZ2j . Here, we investigate stability of periodic traveling waves: specifi- 
cally, sufficient conditions for stability of the wave. Our main result, generalizing results of 
Oh and Zumbrun [QZ4] in dimensions d > 3, is to show that strong spectral stability in the 
sense of Schneider ]S1\ \S2\ implies linearized and nonlinear L 1 n H K — > L°° bounded 
stability, for all dimensions d > 1, and asymptotic stability for dimensions d > 2. 
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More precisely, we show that small L n H s perturbations of a planar periodic solution 
u{x,t) = u(x\) (without loss of generality taken stationary) converge at Gaussian rate in 
L p , P > 2 to a modulation 

(1.1) u(xi-ip(x,t)), 

of the unperturbed wave, where x = (xi,x), x = (x2, ■ ■ ■ , x^), and ip is a scalar function 
whose x- and t- gradients likewise decay at least at Gaussian rate in all L p , p > 2, but which 
itself decays more slowly by a factor i 1 / 2 ; in particular, ip is merely bounded in L°° for 
dimension d = 1. 

The one-dimensional study of spectral stability of spatially periodic traveling waves of 
systems of viscous conservation laws was initiated by Oh and Zumbrun |OZlj in the "quasi- 
Hamiltonian" case that the traveling-wave equation possesses an integral of motion, and in 
the general case by Serre [Sel]. An important contribution of Serre was to point out a larger 
connection between the linearized dispersion relation (the function A(£) relating spectra to 
wave number of the linearized operator about the wave) near zero and the formal Whitham 
averaged system obtained by slow modulation, or WKB, approximation. 

In [OZ3], this was extended to multi-dimensions, relating the linearized dispersion rela- 
tion near zero to 

dtM + Y^dx^i =0, 

(1-2) 

a t (nN) + v x (ns) = o, 

where M G M. n denotes the average over one period, F 3 the average of an associated flux, 
Q = \V X ^\ G M 1 the frequency, S = -*t/|V x *| € M 1 the speed s, and N = V X ^/\S/ X ^\ € 
M. d the normal v associated with nearby periodic waves, with an additional constraint 

(1.3) curl (SIN) = curl = 0. 

As an immediate corollary, similarly as in |OZl| . [Sel| in the one-dimensional case, this 
yielded as a necessary condition for multi-dimensional stability hyperbolicity of the averaged 
system ([L2]) - (fL3]) . 

The present study is informed by but does not directly rely on this observation relating 
Whitham averaging and spectral stability properties. Likewise, the Evans function tech- 
niques used in |Sel|,IOZ3] to establish this connection play no role in our analysis; indeed, the 
Evans function makes no appearance here. Rather, we rely on a direct Bloch-decomposition 
argument in the spirit of Schneider |S1|, IS21 IS3| , combining sharp linearized estimates with 
subtle cancellation in nonlinear source terms arising from the modulated wave approxima- 
tion. The analytical techniques used to realize this program are somewhat different from 
those of [Si t IS2 |, [S3] . however, coming instead from the theory of stability of viscous shock 
fronts through a line of investigation carried out in [OZl[ OZ2"l |QZ3"1 IOZ41 IHoZ] . In partic- 
ular, the nonsmooth dispersion relation at £ = typical for convection-diffusion equations 
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requires different treatment from that of |SH |S2[ [S3] in the reaction diffusion case; see Re- 
mark [2~4l Moreover, we detect nonlinear cancellation in the physical x-t domain rather than 
the frequency domain as in |Sll IS2[ [S3]. The main difference between the present analysis 
and that of [OZ4] is the systematic incorporation of modulation approximation (jl.ip . 



1.1 Equations and assumptions 

Consider a parabolic system of conservation laws 

(1.4) u t + ^2P(u) X] =A x u, 

j 

u € Z//(open) € W 1 , P € M. n , x £ M. d , d > 1, t € R + , and a periodic traveling wave solution 

(1.5) u = u(x ■ v — st), 

of period X, satisfying the traveling-wave ODE u" = (J^ - Vjf 3 (u))' — su' with boundary 
conditions u(0) = u(X) =: uq. Integrating, we obtain a first-order profile equation 

(1.6) u = ^Vjf^u) - su - q, 

j 

where (uo, q, s, i>, X) = constant. Without loss of generality take u = e\, s = 0, so that 
u = u(xi) represents a stationary solution depending only on x\. 

Following |Sel[ IOZ3] IOZ4] . we assume: 

(HI) p G C K+l , K > [d/2] +4. 

(H2) The map H : IxWxIx S^ 1 x W n R n taking (X;a, s,u,q) ^ u(X;a,s,u,q)-a 
is a submersion at point (X; u(Q), 0, e\, q), where u{-; •) is the solution operator for (jl.6p . 

Conditions (H1)-(H2) imply that the set of periodic solutions in the vicinity of u form 
a smooth (n + d + l)-dimensional manifold {u a {x ■ v(a) — a — s(a)t)}, with a £ K, a € R n+d . 



1.1.1 Linearized equations 

Linearizing f j 1 . 4 [) about u(-), we obtain 

(1.7) v t = Lv:=A x v-^2(A j v) X] , 

where coefficients A 3 := Df J (u) are now periodic functions of x\. Taking the Fourier 
transform in the transverse coordinate x = (x2, ■ ■ ■ ,Xd), we obtain 

(1.8) «t = L £V = v xl}X1 - (A 1 v) Xl - i A 3 ijV - ^2 fju, 
where £ = (£2, • • • ,£d) is the transverse frequency vector. 
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1.1.2 Bloch Fourier decomposition and stability conditions 

Following O EH I52j IS3j. we define the family of operators 

(1.9) L ( = < + l.j' iiJ ' : 

operating on the class of I? periodic functions on [0, X]; the (L 2 ) spectrum of is equal 
to the union of the spectra of all with £i real with associated eigenfunctions 

(1.10) w( Xl ,l\) :=e* ia,1 9(si,£i,£A), 

where q, periodic, is an eigenfunction of Lg. By continuity of spectrum, and discreteness of 
the spectrum of the elliptic operators on the compact domain [0, X], we have that the 
spectra of may be described as the union of countably many continuous surfaces Aj(£). 
Without loss of generality taking X = 1, recall now the Bloch-Fourier representation 

(±Y f" I 



(1.11) „( l) = (_j / / e^'MS, Xl)d(i d( 



-1 



of an I? function u, where u(£,xi) := ^2 k e 2mkxi u(^i + 27r/c, £) are periodic functions of 
period X = 1, u(£) denoting with slight abuse of notation the Fourier transform of u in the 
full variable x. By Parseval's identity, the Bloch-Fourier transform u{x) — > u(£,xi) is an 
isometry in L 2 : 

(1-12) IMI-L 2 ^) = ||u|| L 2( ?;i 2( ;Cl )), 

where L 2 (x\) is taken on [0,1] and L 2 {^) on [— 7r,7r] x Moreover, it diagonalizes the 

periodic-coefficient operator L, yielding the inverse Bloch-Fourier transform representation 

(1.13) e Lt u = (1-Y [ W [ e^e^MCx^di 

relating behavior of the linearized system to that of the diagonal operators L^. 

Following [OZ4j . we assume along with (H1)-(H2) the strong spectral stability conditions: 
(Dl) a(L£ C {ReA < 0} for £ ^ 0. 

(D2) Rea(L ? ) < -#|£| 2 , 9 > 0, for £ G R d and |f | sufficiently small. 
(D3) A = is a semisimple eigenvalue of Lq of multiplicity exactly n + 10 
For each fixed angle £ := £/|£|, expand = Lq + l^l-L 1 + |£| 2 £ 2 - By assumption (D3) 
and standard spectral perturbation theory, there exist n + 1 smooth eigenvalues 

(1.14) A i (0 = -m i (0 + Kiei) 



1 The zero eigenspace of Lo is at least (n + l)-dimensional by the linearized existence theory and (H2), 
and hence n + 1 is the minimal multiplicity; see Scl ( OZ3]. As noted in [OZ11 [OZ3] . minimal dimension 
of this zero eigenspace implies that (M, NQ) of (|1.2I) gives a nonsingular coordinatization of the family of 
periodic traveling-wave solutions near u. 
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of bifurcating from A = at £ = 0, where — iaj are homogeneous degree one functions 
given by |£| times the eigenvalues of Tl$L l lKerL ' w ^ n ^0 * ne zero eigenprojection of Lq. 

Conditions (D1)-(D3) are exactly the spectral assumptions of [S11IS2| I53]. corresponding 
to "dissipativity" of the large-time behavior of the linearized system. As in [OZ4] . we make 
the further nondegeneracy hypothesis: 

(H3) The eigenvalues A = — icij(£)/|£| of iloL^- erL are simple. 
The functions a,- may be seen to be the characteristics associated with the Whitham av- 
eraged system (|1.2p - (|1.3jl linearized about the values of M, S, N, 0, associated with the 
background wave u; see [OZ3| [QZ4] . Thus, (Dl) implies weak hyperbolicity of (|1.2p - (|1.3|) 
(reality of a,), while (HI) corresponds to strict hyperbolicity. 

1.2 Main results 

With these preliminaries, we can now state our main results. 

Theorem 1.1. Assuming (Hl)-(HS) and (Dl)-(DS), for some C > and ip G W K '°°(x,t), 

\u - «(• - i>)\ L p{t) < c(i + tyi^-yp)\u - u\ LlnH K\ t =o, 

(1.15) \u-u(- - ip)\ H K{t) < C(l + t)~*|u- u\ L i nH K\ t =o, 

\(ik,il> x )\ W K + ia, < C(l + ty^ l - l M\u - u\ L i nH K\ t=0 , 

and 

(1.16) \u-u\ LP {t), \^{t)\ L p < C(\ + t)-i {l ~p )+ ^\u-u\ L i nHK \t=o 

for all t > 0, p > 2, d = 1, for solutions u of (jl.4p with \u — u\ L i nH K\ t= o sufficiently small. 
In particular, u is nonlinearly bounded L 1 n H K — > L°° stable for dimension d = 1. 

Theorem 1.2. Assuming (Hl)-(HS) and (Dl)-(DS), for any e > 0, some C > and 

|u - < - ^)\ L v{t) < C(l + tj-ff 1 - 1 /?)^ - n| z i nH K| t=0 , 

(1.17) |u-u(- -^)lij-if(t) < C^ + ^'^lw-uliin^K^o, 

|(Vt, ^)|^+i,p < C(l + t)-fCi-i/i»)+=-i| fi _ u| z i nH K | t=0 , 

and 

|u-u| LP (t), I^WIlp <C(l + t)-^ (1 -p )+£ |n-n| L i nH K|i=o, 
|tt-u| H K(t), I^WIirx < C(l + t) i +e \u - u\ L i nHK \ t =o, 

for all t > 0, p > 2, d = 2, for solutions u of (jl.4p with \u — u\ L i nH K\t=o sufficiently small. 
In particular, u is nonlinearly asymptotically L 1 D H K — > H stable for dimension d = 2. 
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Remark 1.1. In Theorem derivatives in x £ M 2 refer to total derivatives. Moreover, 
unless specified by an appropriate index, throughout this paper derivatives in spatial variable 
x will always refer to the total derivative of the function. 

In dimension one, Theorem 11.11 asserts only bounded L 1 n H — > L°° stability, a very 
weak notion of stability. The absence of decay in perturbation u — u indicates the delicacy 
of the nonlinear analysis in this case. In particular, it is crucial to separate off the slower- 
decaying modulated behavior (jl.ip in order to close the nonlinear iteration argument. 

Remark 1.2. In dimension d = 1, it is straightforward to show that the results of Theorem 
11.11 extend to all 1 < p < oo using the pointwise techniques of [ OZ2| ; see Remark 13.31 

Remark 1.3. The slow decay of \u — u|j>(£) ~ \4>{t)\Lp in f|l. 16[) is due to nonlinear 
interactions; as shown in |OZ2tlOZ4] . the linearized decay rate is faster by factor (l + t) -1 / 2 
(Proposition 12.1] ). In [QZ4] . it was shown that for d > 3, where linear effects dominate 
behavior, (|1.16p may be replaced by the stronger estimate \u — u|lp(£), \ip(t)\LP < C(l + 

t) 2 p'\u — u\ L i nH K\t=o- These distinctions reflect fine details of both linearized estimates 
(Section [3|) and nonlinear structure (Sections l4.lH4.2p that are not immediately apparent 
from the formal Whitham approximation ()1.2[) — ()1.3[) . 

I. 3 Discussion and open problems 

Linearized stability under the same assumptions, with sharp rates of decay, was established 
for d = 1 |OZ2] and for d > 1 in [OZ4J, along with nonlinear stability for d > 3. Theorem 

II. 11 completes this line of investigation by establishing nonlinear stability in the critical 
dimensions d = 1, 2, a fundamental open problem cited in |OZH [OZ4J. 

This gives a generalization of the work of |Sll IS2j. IS3| for reaction diffusion equations 
to the case of viscous conservation laws. Recall that the analysis of (S 1 1. |S2} [S3] concerns 
also multiply periodic waves, i.e., waves that are either periodic or else constant in each 
coordinate direction. It is straightforward to verify that the methods of this paper apply 
essentially unchanged to this case, to give a corresponding stability result under the analog 
of (H1)-(H3), (D1)-(D3), as we intend to report further in a future work. Likewise, the 
extension from the semilinear parabolic case treated here to the general quasilinear case is 
straightforward, following the treatment of [OZ4]. 

On the other hand, as noted in [OZ2j . condition (D3) is in the conservation law setting 
nongeneric, corresponding to the special "quasi-Hamiltonian" situation studied there; in 
particular, it implies that speed is to first order constant among the family of spatially pe- 
riodic traveling-wave solutions nearby u. In the generic case that (D3) is violated, behavior 
is essentially different |OZH [QZ2] . and perturbations decay more slowly at the linearized 
level. Nonlinear stability remains an interesting open problem in this setting. 

Our approach to stability in the critical dimensions d = 1, 2, as sugg ested in [OZ4 , is, 
loosely following the approach of |Sll IS2l IS3] , to subtract out a slower-decaying part of the 
solution described by an appropriate modulation equation and show that the residual decays 
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sufficiently rapidly to close a nonlinear iteration. It is worth noting that the modulated 
approximation u{x\ — ip(x, t)) of f 1 1 . 1 j) is not the full Ansatz 



^(x, t) := xi — ifi(x, t), associated with the Whitham averaged system (|1.2p - (|1.3p . where u a 
is the manifold of periodic solutions near u introduced below (H2), but only the translational 
part not involving perturbations a in the profile. (See [OZ3 for the derivation of (|1 .2[) 
(jl.3p and (|1.19p .) That is, we don't need to separate out all variations along the manifold 
of periodic solutions, but only the special variations connected with translation invariance. 

The technical reason is an asymmetry in y-derivative estimates in the parts of the Green 
function associated with these various modes, something that is not apparent without a 
detailed study of linearized behavior as carried out here. This also makes sense formally, 
if one considers that (|1.2p indicates that variables a, V a . v I / are roughly comparable, which 
would suggest, by the diffusive behavior ^ >> V x ^, that a is neglible with respect to 

However, note that in the case that (D3) holds, hence wave speed is stationary along the 
manifold of periodic solutions, the final equation of (ll.2p decouples to (^a;)* = (&N)t = 0, 
and could be written as = in terms of \I/ alone. Hence, there is some ambiguity in this 
degenerate case which of ^, fy x is the primary variable, and in terms of linear behavior, the 
decay of variations a and ^ are in fact comparable [OZ4]; in the generic case, a and ^f x are 
comparable at the linearized level [OZ2] . It would be very interesting to better understand 
the connection between the Whitham averaged system (or suitable higher-order correction) 
and behavior at the nonlinear level, as explored at the linear level in |OZ3t |QZ4"1 \JZ1\ DZB . 

2 Basic linearized stability estimates 

We begin by recalling the basic linearized stability estimates derived in [OZ4] , We will 
sharpen these afterward in Section [3l By standard spectral perturbation theory [K], the 
total eigenprojection P(£) onto the eigenspace of L^ associated with the eigenvalues Aj(£), 
j = 1, . . . , described in the introduction is well-defined and analytic in £ for £ sufficiently 
small, since these (by discreteness of the spectra of L^) are separated at £ = from the rest 
of the spectrum of Lq. Introducing a smooth cutoff function that is identically one 
for |£| < e and identically zero for |£| > 2e, e > sufficiently small, we split the solution 
operator S(t) := e Lt into low- and high-frequency parts 



(1.19) 



u a (^(x,t)) 



(2.1) 




and 



(2.2) 
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2.1 High-frequency bounds 

By standard sectorial bounds |Hel IPaj and spectral separation of Aj(£) from the remaining 
spectra of Lg, we have trivially the exponential decay bounds 

|!eV(/_^))/|| LW]) < Ce- et \\f\\ LH[0tX]) , 
(2-3) \\e L ^(I - (f>P(0)d l x J\\ L 2 {[0tX]) < Ct-h- dt \\f\\ Lmx]) , 

\\d l xl e L t\l- cf>P(0)f \\ LH[0 ,x]) < Ct- L ie- et \\f\\ LH[0 , x]) , 

for 9, C > 0, and < m < K (K as in (HI)). Together with (11.121) . these give immediately 
the following estimates. 

Proposition 2.1 QOZ4]). Under assumptions (H1)-(H3) and (D1)-(D2), for some 9, 
C> 0, and all t > 0, 2 < p < oo, < / < K + 1, < m < K , 

\\d x s u (t)f\\ mx) , \\s n {t)$j\\v {x) < ct-h- et \\f\\ LHx) , 

\\d?S H {t)f\\ LP{x) , ||S IJ (i)<9£7|| iP(x) < Ct-^-^e- dt \\f\\ L2{x) , 

where, again, derivatives in the variable x € M. d refer to total derivatives. 

Proof. The first inequalities follow immediately by (|1.12p and (|2.3p . The second follows for 
x\ derivatives in the case p = oo, m = by Sobolev embedding from 

\\s n (t)f\\ Lao[i . LHxi)) < cr^e- et \\f\\ L2{[0>x]) 

and 

Wdx^ 11 {t)f\\L°°(x;L2(xi)} < Ct 4 2e _6 '*||/|| L 2([ 0iX ]), 

which follow by an application of (|1.12p in the x\ variable and the Hausdorff- Young in- 
equality < ||/|lx,x(|) i n the variable x. The result for derivatives in x\ and general 
2 < p < oo then follows by L p interpolation. Finally, the result for derivatives in x follows 
from the inverse Fourier transform, equation (|2.2p . and the large |£| bound 

\e U f\L*(xi) < z~ 6 ^ H \f\l?(x x ), Ifl sufficiently large, 

which easily follows from Parseval and the fact that Lg is a relatively compact perturbation 
of d x — |£| 2 . Thus, by the above estimate we have 

\\e Lt d x f\\ L 2 (x) < C\\e L ^\i\f\\ L 2 {xi>0 

<Csup (e-^l 2 *|C|) 11/11^(^,0 

A similar argument applies for 1 < m < K. □ 
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2.2 Low-frequency bounds 

Denote by 

(2.5) G I (x,t-y):=S I (t)8 y (x) 
the Green kernel associated with S , and 

(2.6) [G|(xi,t;yi)] := 0(O^(Oe Let [^(xi)] 

the corresponding kernel appearing within the B loch-Fourier representation of G 1 , where 
the brackets on [Gt] and [S y ] denote the periodic extensions of these functions onto the 
whole line. Then, we have the following descriptions of G 1 , [G|], deriving from the spectral 
expansion (|1.14p of Lg near £ = 0. 

Proposition 2.2 ( [QZ4] ). Under assumptions (Hl)-(HS) and (D1)-(D3), 

n+1 

[Gf(x 1 ,t;y 1 )]=0(OX)e A itt) t ?J -(^x 1 )^(e,yi)* 1 

3=1 

(2.7) G\x,t ]y ) = {^) A J S< x - y \Gl{ Xl ,t- yi )n 

where * denotes matrix adjoint, or complex conjugate transpose, </j(£, •) and gj(£, •) are right 
and left eigenf unctions of L^ associated with eigenvalues Aj(£) defined in (|1.14p . normalized 
so that (q~j,qj) = 1, where Aj/|£| is a smooth function of |£| and £ := £/|£| and g^ and gj 
are smooth functions of |£|, £ := and or yi, imt/j 5RAj(£) < — #|£| 2 . 

Proof. Smooth dependence of Xj and of g, q as functions in L 2 [0,X] follow from standard 
spectral perturbation theory [K] using the fact that Xj split to first order in |£| as £ is varied 
along rays through the origin, and that L^ varies smoothly with angle £. Smoothness of 
qj, q~j in x±, y\ then follow from the fact that they satisfy the eigenvalue equation for L^, 
which has smooth, periodic coefficients. Likewise, (|2.7|) (i) is immediate from the spectral 
decomposition of elliptic operators on finite domains. Substituting (|2.5j) into (|2.1h and 
computing 

(2.8) 8y(£,xi) =J2 e2 ™ kxi $y(€ + 27Tke i) = e 2 * ikxi e-* y - 2 * ikyi = e^^^i)], 

k k 

where the second and third equalities follow from the fact that the Fourier transform either 
continuous or discrete of the delta-function is unity, we obtain 

G I (x,t;y) = (±-Y f / e^PiQ^SyfoxJdt 

V27T/ J-TrjRd-1 
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yielding (|2.7p (ii) by (|2.6p (i) and the fact that 4> is supported on [— tt, it]. □ 
Proposition 2.3 ( [OZ4] ). Under assumptions (Hl)-(HS) and (Dl)-(DS), 

(2.9) sup||Gr J (-,t,;j,)[| ii p (as)) sup R )2/ G 7 (-, t, ; y)\\ LP(x) < G(l + ty^S 

y y 

for all 2 < p < oo, t > 0, where C > is independent of p. 

Proof. Prom representation (12.7p (ii) and 3iAj(£) < — 0\£,\ 2 , we obtain by the triangle in- 
equality 

(2.10) II^IUoc^) < c||e- e ^ a V(Olki(o < c(i + t)-i, 

verifying the bounds for p = oo. Derivative bounds follow similarly, since derivatives falling 
on qj or qj are harmless, whereas derivatives falling on e^'^^' bring down a factor of £, 
again harmless because of the cutoff function <j>. 

To obtain bounds for p = 2, we note that (|2.7p (ii) may be viewed itself as a Bloch- 
Fourier decomposition with respect to variable z := x — y, with y appearing as a parameter. 
Recalling (I1.12p . we may thus estimate 

sup 1 1 G J (a?,*;y) | \ L i {x) = V" sup W^e^^q^-, z{)^(; yi)h^^(zie[o,X])) 
y j y 

( 2 - n ) < ^y]sup||(/)(0e^ l5|2 '|| L 2 (c) ||g i || L 2 (0iX) ||q ; i || L o O(0iX) 

3 V 
< C(l + t)~i, 

where we have used in a crucial way the boundedness of g,-; derivative bounds follow simi- 
larly. Finally, bounds for 2 < p < oo follow by L p -interpolation. □ 

Remark 2.4. In obtaining the key L 2 -estimate, we have used in an essential way the 
periodic structure of qj, qj. For, viewing G 1 as a general pseudodifferential expression 
rather than a Bloch-Fourier decomposition, we find that the smoothness of qj, <jj is not 
sufficient to apply standard L 2 — > I? bounds of Hdrmander, which require blowup in £ 
derivatives at less than the critical rate found here; see, e.g., [Hj for further discussion. 
Nor do the weighted energy estimate techniques used in [Sl[|S2"l lS3] apply here, as these also 
rely on the property of smoothness of Xj, qj, qj with respect to £ at the origin £ = 0. The 
lack of smoothness of the linearized dispersion relation at the origin is an essential technical 
difference separating the conservation law from the reaction diffusion case; see jOZ4j for 
further discussion. 

Remark 2.5. Underlying the above analysis, and also the technically rather different 
approach of [OZ2], is the fundamental relation 

(2.12) G(x,t;y) = (±-Y f / e^-^G^i^M 

V27T/ J-TrjRd-l 
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which, provided cr(L^) is semisimple, yields the simple formula 

G(x,t-y) = (±) d f / e^-y^^ ]t ^^)^yiYdi 

resembling that of the constant-coefficient case, where Xj runs through the spectrum of L^. 
The basic idea in both cases is to separate off the principal part of the series involving small 
Aj(£) and estimate the remainder as a faster-decaying residual. 

Corollary 2.6 f [OZ4] ). Under assumptions (Hl)-(HS) and (Dl)-(DS), for allp >2,t>0, 

(2.i3) hs'w/Hlp, ||9x5 7 (t)/|| LP , iis^r/Hw < c(i + tr* (1 -%iUx. 

Proof. Immediate, from (|2.9p and the triangle inequality, as, for example, 

\\S I (t)f(-)\\ I ,= / G I (x,t;y)f(y)dy _ < / supWG 1 (;t;y)\\ LP \f(y)\dy. 



Lp{x) 



U 



□ 



Proposition 2.1 ([QZ4]). Assuming (Hl)-(HS), (Dl)-(DS), for some C > 0, all t > ; 

P > 2, < I < K, 

(2.14) \\S(t)d l x u \\LP <cH(l + t)"^l4)+^-f-||| U0 || LlnIj2 . 

Proof. Immediate, from (|2.4p and (|2.13p . □ 

2.3 Additional estimates 

Lemma 2.7. Assuming (H1)-(H3), (D1)-(D3), for all t > 0, < I < K , 

(2.15) w^mi^, ||5 7 (t)^/|| iP(x) ^cii + tyi^-^wfw^. 

Proof. Prom boundedness of the spectral projections Pj(£,) = Qj(Qj, •) in £ 2 [0,^] and their 
derivatives, another consequence of first-order splitting of eigenvalues Aj(£) at the origin, 
we obtain boundedness of </>(£)-P(£) eLf;t an d thus, by (|1.12p . the global bounds 

(2.16) wdi&miL^), ws'^difu^ < c\\f\\ LHx) , 

for all t > 0, yielding the result for p = 2. Moreover, by boundedness of q, q in all L p (xi), 
we have 

mnOe Lst f&-)\L<*>( Xl ) < Ce-W H \P(Of(ti,-)\L°°( xl) < Ce- e ^ 2t \f(C,-)\ LHxi) , 
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C, 6 > 0, yielding by S 1 f = / R a-i e^(£)P(Oe^/(£^iRi df the bound 

(2.17) * (hYLL-, cme^iH^^di 

<c\m^ m2 \^)\f\L^, xl) 

= C(l + t)-i\\f\\ L 2 {[0iX]) , 

yielding the result for p = oo, Z = 0. The result for p = oo, 1 < Z < K follows by a similar 
argument. The result for general 2 < p < oo then follows by L p interpolation between p = 2 
and p = oo. □ 

By Riesz-Thorin interpolation between (|2.15j) and (12. 13ft . we obtain the following, ap- 
parently sharp bounds between various L q and LP @ 

Corollary 2.8. Assuming (HO)-(HS) and (Dl)-(DS), for all 1 < q < 2 < p, t > 0, 

0<1<K, 

(2.18) ||^(t)/|| iP , ||^(t)ai/|| LP <C(l + t)-^ ( ^^||/|| L9 . 

Proposition 2.2. Assuming (Hl)-(HS), (Dl)-(DS), for some C > 0, all t > 0, 1 < q < 

2<p, and < I < K, 

(2.19) \\s(t)dluo\\ LP < C(i + t)--^-p + n--^-* ) ->\\ 

Proof. Immediate, from (|2,4p and (|2.8p , □ 



3 Refined linearized estimates 

The bounds of Proposition 12.11 are sufficient to establish nonlinear stability and asymptotic 
behavior in dimensions d > 3, as shown in [OZ4]. However, they are not sufficient in the 
critical dimensions d = 1,2; see Remark 1, Section 7 of [OZ4 . Comparison with standard 
diffusive stability arguments as in |Z7] show that this is due to the fact that the full solution 
operator \S(t)d x \ decays no faster than S(t), or, equivalently, G y no faster than G. 

Following the basic strategy introduced in |ZFH \Z1\ IMaZ2} IMa Z4] in the context of vis- 
cous shock waves, we now perform a refined linearized estimate separating slower-decaying 
translational modes from a faster-decaying "good" part of the solution operator. This will 
be used in Section |4] in combination with certain nonlinear cancellation estimates to show 
convergence to the modulated approximation (jl.ip at a faster rate sufficient to close the 
nonlinear iteration. 

The key to this decomposition is the following observation. 

2 The inclusion of general p > 2 in Lemma [2.71 repairs an omission in |OZ4] . where the bounds (|2.8|l were 
stated but not used. 
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Lemma 3.1. Assuming (Hl)-(HS), (Dl)-(DS), let A,-(£/|£U)> ?;(£/|£U,-), &(£/|£l,& 
denote the eigenvalues and associated right and left eigenf unctions of L^, with qj, <jj smooth 
functions o/£/|£| and |£| as noted in Prop. \2.S\ . Then, without loss of generality, qi(ui, 0, •) = 
v! , while q~j(u), 0, •) for j ^ 1 are constant functions depending only on angle uj = 

Proof. Expanding = Lq + l^l-^/^i + |£| 2 L 2 ^| as in the introduction, consider the con- 
tinuous family of spectral perturbation problems in |£| indexed by angle u = Then, 
both facts follow by standard perturbation theory [K] using the observations that v! is in 
the right kernel of Lq and constant functions c are in the left kernel of Lq, with 

(c,L l v!) = (c, (ujx{2d Xl - A x ) - ^u)jAj))u') = (c^xd^u -^u)jd Xl f J (u)) = 0, 

where (•,•) denotes L 2 (x\) inner product on the interval x\ € [0, X], that the dimension 
of kerLo by assumption is (n + 1), so that the orthogonal complement of v! in KerLo 
is dimension n so exactly the set of constant functions, and that by (H3) the functions 
qj(u, 0, •) and qj(ui, 0) are right and left eigenfunctions of IIoL 1 |kerLo (Ho as earlier denoting 
the zero eigenprojection associated with Lq). □ 

Remark 3.2. The key observation of Lemma 13.11 can be motivated by the form of the 
Whitham averaged system (|1 .2j) . For, recalling (Section 1 1 . 3|) that (D3) implies that speed 
s is stationary to first order at u along the manifold of nearby periodic solutions, we find 
that the last equation of (|1.2p reduces to (Vjf )t = 0, i.e., the equation for the translational 
variation ^ decouples from the equations for variations in other modes. This corresponds 
heuristically to the fact derived above that the translational mode u'{x\) decouples in the 
first-order eigenfunction expansion. 

Corollary 3.1. Under assumptions (H1)-(H3), (Dl)-(DS), the Green function G(x,t;y) 
of (]1 . 7H decomposes as G = E + G, 

(3.1) E = u'(x)e(x,t;y), 

where, for some C > 0, all t > 0, i<q<2<p< oo, < j, k,l, j + 1 < K , 1 < r < 2, 

/+oo 
G{x,t-y)f{y)dy < C(l + t)"f (i/a-i/rt t - 5 d/«-i/2) |/| L q ^ 
-oo LP (°>) 



(3.2) 



-oo 
rOO 

d r y G(x,t;y)f(y)dy < C(l + t)-*<V2-i/iO- 3 + 5 

xr l(l/ 9 -l/ 2 )^ l/IL9ni2! 

hoo 

%G(x,t;v)f(y)dy < C(l + t )- 5(1/2-1/*)- 5 +r 



xdw^-i/i 



L«nL 2 • 
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(3.3) / did^e(x,t;y)f(y)dy < (1 + tJ-sd/ff-Vrt-V 1 |/| i9 . 



Moreover, e(x,t;y) = /or t < 1. 

Proof. We first treat the simpler case 9 = 1. Recalling that 

(3.4) G 7 (x,t;y) = (-) / e^UiO £ e A ^(£, 

VZ7r/ JR d . =1 

define 

(3.5) e(x,t;y) = (^-V / e^-^(£) e MOt^ yi)X) 

so that 
(3.6) 



n+l 



G'foijy) -u'(x 1 )e(x,t;y) = (±-Y [ e*<*^U(o£^\(Z/\t\Axi)m,Vi)*<% 

+ (2^) J d E e^-^)0(e)e A ^«*O(|e|)^- 
Noting, by Lemma 13. 11 that d y q(uj,0,y) = constant for j ^ 1, we have therefore 

(3.7) d r y (G I (x,t;y)-u , (x 1 )e(x,t;y)) = (— ) jf^ e*(*-»V(0 ^ e A '«>'0(|£|K, 
which readily gives 

(3.8) |^(G 7 (x, i; y) - u'( Xl )e(x, t; y))\ Lv < C(l + t)"f (^/pH, 
p > 2, by the same argument used to prove (|2.9|) . and similarly 

(3.9) \dl{G\x, t; y) - u'(x 1 )e(x, t; y))\ LP < c(l + t)~i (^K. 



These yield ()3.2p by the triangle inequality. 

Defining e(x, t; y) := x(*)e(x, y)> where x is a smooth cutoff function such that x(t) = 1 
for t > 2 and x(i) = for i < 1, and setting G := G — u'(xi)e(x, t; y), we readily obtain the 
estimates flO} by combining fl3ji with bound ([23]) on G 7/ . Bounds fl33D follow from P^fl 
by the argument used to prove (|2.9p . together with the observation that x- or t-derivatives 
bring down factors of followed again by an application of the triangle inequality. 

The cases 1 < q < 2 follow similarly, by the arguments used to prove (|2.15p and (|2.8p . □ 
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Remark 3.3. Despite their apparent complexity, the above bounds may be recognized 
as essentially just the standard diffusive bounds satisfied for the heat equation [Z7J. For 
dimension d = 1, it may be shown using pointwise techniques as in [OZ2] that the bounds 
of Corollary 13.11 extend to all 1 < q < p < oo. 

Note the strong analogy between the Green function decomposition of Corollary 13,11 
and that of [MaZ3j IZ4] in the viscous shock case. We pursue this analogy further in the 
nonlinear analysis of the following sections, combining the "instantaneous tracking" strategy 
of |ZHl IZll IZ41 |Z7, MaZ2, MaZ4] with a type of cancellation estimate introduced in [HoZj. 

4 Nonlinear stability in dimension one 

For clarity, we carry out the nonlinear stability analysis in detail in the most difficult, 
one-dimensional, case, indicating afterward by a few brief remarks the extension to d = 2. 
Hereafter, take x € Mr, dropping the indices on p and Xj and writing u t + f(u) x = u xx . 

4.1 Nonlinear perturbation equations 

Given a solution u(x,t) of (|1.4p . define the nonlinear perturbation variable 

(4.1) v = u — u = u(x + tp(x, t)) — u(x), 
where 

(4.2) u(x,t) := u(x + tp(x,t)) 
and ip : R x R — >• R is to be chosen later. 

Lemma 4.1. For v, u as in (j4.ip . (|4.2j) . 

(4.3) u t + f{u) x - u xx = (d t - L) u'{ Xl )^{x, t) + 8 X R + (d t + d 2 x )S, 
where 

R := win + v^ xx + (u x + %)t^7T = °(l u l(l^tl + I^D + (^fZ^Tf ) 

and 

S ■= -vip x = 0(\v\(\ip x \). 

Proof. To begin, notice from the definition of u in (|4.2|) we have by a straightforward 
computation 

u t (x,t) = u x (x + ip(x,t),t)ijj t (x,t) + u t (x + 1p,t) 
f(u(x,t)) x = df(u(x + ip(x,t),t))u x (x + ip,t)-(l + if> x (x,t)) 
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and 

u xx (x,t) = (u x (x + ip(x,t),t) ■ (1 + ip x (x,t))) x 

= u xx (x + ip(x,t),t) ■ (l + ip x (x,t)) + (u x {x + if>(x,t),t) ■ ip x (x,t)) x . 

Using the fact that ut + df(u)u x — u xx = 0, it follows that 

^ ^ u t + f(u) x - u xx = u x ip t + df(u)u x ip x - u xx ip x - {u x ip x ) x 

= u x ip t - U t lf) x - (u x tp x ) x 

where it is understood that derivatives of u appearing on the righthand side are evaluated 
at (x + ip(x,t),t). Moreover, by another direct calculation, using the fact that L(u'(x)) = 
by translation invariance, we have 

(d t - L) u'(x)ip = u x ip t ~ utipx ~ {u x ij) x ) x . 

Subtracting, and using the facts that, by differentiation of (u + v)(x, t) = u(x + ip,t), 

U x + V x = U x (l + 7p x ), 

u t + v t = u t + u x ip t , 



(4.5) 
so that 



/ — i s 

u x -u x -v x = -{u x + V x ) 



(4.6) 



we obtain 



u t -u t -v t = -(u x + v x ) 



l + lpx' 

A 

1 + ^' 



Ut + f(u) x - u xx = (d t - L)u'(x)ip + v x i>t - v t ipx - (v x il>x)x + {{u x + v x ) ^ x ) , 

V 1 + 1p x J X 

yielding (g^J by v x ip t ~ v t ip x = {wip t )x ~ (vip x )t and (v x ip x ) x = (vip x ) xx - (vip xx ) x . □ 
Corollary 4.2. The nonlinear residual v defined in (|4.ip satisfies 

(4.7) v t -Lv = (d t - L) u\x x )i> -Q x + R x + (d t + 8 2 X )S, 
where 

(4.8) Q := f(u(x + 4,(x, t),t)) - f{u{x)) - df(u(x))v = 0(\v\ 2 ), 

4> 2 

(4.9) R := vip t + vipxx + (u x + v x ) \* , 
and 

(4.10) S :=-viP x = 0(\v\(\^ x \). 

Proof. Taylor expansion comparing (|4.3|) and u\ + f(u) x — u xx = 0. □ 
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4.2 Cancellation estimate 

Our strategy in writing (|4.7|) is motivated by the following basic cancellation principle. 
Proposition 4.3 ( |HoZ] ). For any f(y,s) G IP n C 2 with f(y,0) = 0, there holds 

(4.11) J* J G{x, t - s; y)(d s - L y )f(y, s)dy ds = f{x, t). 

Proof. Integrating the left hand side by parts, we obtain 
(4.12) 

G{x, 0; y)f(y, t)dy - [ G{x, t; y)f(y, 0)dy + f I {d t - L y )*G{x, t - s; y)f{y, s)dy ds. 



Noting that, by duality, 

{d t - L y )*G(x, t- s;y) = 6{x - y)5(t - s), 

5{-) here denoting the Dirac delta-distribution, we find that the third term on the righthand 
side vanishes in (|4.12p . while, because G(x, 0; y) = 5{x — y), the first term is simply f{x,t). 
The second term vanishes by f{y, 0) = 0. □ 

Remark 4.1. For ip = tp{t), term {dt — L)u'ip in (|4.7|) reduces to the term ip{t)u'{x) 
appearing in the shock wave case [ZHl [ZTJ EU EZ1 IMaZ2l IMaZ4] , 

4.3 Nonlinear damping estimate 

Proposition 4.2. Let vq G H k (K as in (HI)), and suppose that for < t < T, the H K 

norm of v and the H K {x,t) norms of ipt an d remain bounded by a sufficiently small 
constant. There are then constants 9\2 > so that, for all < t <T, 

(4.13) \v{t)\ 2 HK < Ce- e ^\v{0)\ 2 HK + C j\- e ^ (|«|| 2 + |(^, Vfe)lfr* (a ,,t)) W da. 

Proof. Subtracting from the equation (|4.4p for u the equation for u, we may write the 
nonlinear perturbation equation as 

(4.14) v t + {df{u)v) x - v xx = Q{v) x + u x ip t - u t ip x - {u x ip x ) x , 

where it is understood that derivatives of u appearing on the righthand side are evaluated 
at {x + ip{x, t),t). Using (|4.6p to replace u x and ut respectively by u x + v x — {u x + v x ) 1 ^ 

and lit + vt — {u x + v x ) 1 ^ , and moving the resulting vtip x term to the lefthand side of 
([4TIi]) . we obtain 

(1 + ip x )v t - v xx = -{df{u)v) x + Q{v) x + u x ip t 

(4-15) / V 2 \ 

- {{u x + v x )ip x ) x + [{u x +v x )— — - . 

V 1 + tp x / X 
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Taking the L 2 inner product in x of Ylf=o against (|4.15p . integrating by parts, and 
rearranging the resulting terms, we arrive at the inequality 

dt\v\ 2 HK {t) < -e\d« +1 v\ 2 L2 + c (\v\ 2 HK + mM\ 2 H K( X ,t)) > 

for some 6 > 0, C > 0, so long as \u\ h k remains bounded, and \v\ h k and \(tptj'4 , x)\H K (x,t) 
remain sufficiently small. Using the Sobolev interpolation \v\ 2 HK < Id^^v]^ + C\v\ 2 L2 for 
C > sufficiently large, we obtain d t \v\ 2 HK {t) < -6\v\ 2 HK + C (\v\\ 2 + \(t/j t , i>x)\ 2 H K( x ^ 
from which (14. 131) follows by Gronwall's inequality. □ 

4.4 Integral representation/ ^-evolution scheme 

By Proposition H31 we have, applying Duhamel's principle to (|4.7|) . 

poo 

v(x,t) = I G(x,t;y)v (y)dy 

t roo 



(4.16) ■ x 



+ 



ft f'OO 

/ / G(x,t-s;y)(-Qy + R x + St + Syy)(y,s)dyds + ip(t)u'(x). 

JO J-oo 



Defining ip implicitly as 
(4.17) 



oo 

ip(x,t) = -l e(x,t;y)u (y)dy 

oo 

t r+oo 



ft p+oo 

/ / e(x,t-s;y)(-Q y + R x + S t + Syy)(y,s)dyds, 

JO J-OD 



following [ZHl [H IMaZ2l IMaZ3j . where e is defined as in ([53]) . and substituting in (|4~TB|> 
the decomposition G = u'(x)e + G of Corollary 13.11 we obtain the integral representation 



oo 



v(x,t) = I G(x,t;y)v (y)dy 
t roo 



(4.18) 



ft roo 

+ G(x,t- s;y)(-Q y + R x + S t + S yy )(y,s)dyds, 

Jo J -oo 

and, differentiating (|4. 17[) with respect to t, and recalling that e(x,s;y) = for s < 1, 



&> t d^{x,t) = - I d 1 t d^e(x,t;y)u (y)dy 

(4.19) 



oo 

ft f + OO 

/ / d{d^e(x,t- s;y)(-Q y + R x + S t + S yy )(y,s)dyds. 

Jo J-oo 



Equations (|4.18p . (|4.19j) together form a complete system in the variables (v, dfip, d x ip), 
0<j<l,0<k<K, from the solution of which we may afterward recover the shift ip via 
(|4.17p . From the original differential equation (|4.7p together with (|4.19p . we readily obtain 
short-time existence and continuity with respect to t of solutions (v,tpt,tp x ) G H K by a 
standard contraction- mapping argument based on (|4.13j) . (|4.17j) . and and (|3.3p . 
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4.5 Nonlinear iteration 

Associated with the solution (u,^ t ,ip x ) of integral system (|4,18|) - (|4.19|) . define 
(4 20) ({*)'■= SU P \{ v ^t^x)\HK{s){\ + s) 1/4 . 

Lemma 4.3. For all t > for which ((t) is finite, some C > 0, and Eq := \uo\ L i nH K , 

(4.21) at)<c(E +at) 2 ). 

Proof. By (l49j) - (l4T0ll and definition (14201) . 

(4.22) \(Q,R,S)\ L i nL oo < \(v, v x , fa, + \(v, v x , ifr, ^ x )\ 2 Laa < C((t) 2 (l + t)~K 

so long as \if) x \ < \ip x \ H K < remains small, and likewise (using the equation to bound t 
derivatives in terms of x-derivatives of up to two orders) 

(4.23) \(d t + dl)S\ L i nL oo < \(v^ x )\ 2 H2 + 1(^,^)1^,00 <CC(t) 2 (l + t)-i 

Applying Corollary 13.11 with q = 1, d = 1 to representations (|4.18p - ([4.19p . we obtain for 
any 2 < p < oo 

H-,t)\ LPix) <c(i+t)-^-yp) Eo 

+ c((t) 2 [ (i + t-sy^-^it-sy-^i + syUs 

Jo 

icciEo+atf^i+t)-^ 1 - 1 ^ 



(4.24) 



and 
(4.25) 

\{thM(;t)\ WKtP < c{\ + t)-7Eo + cat) 2 [ (i + 1 - s )-K 1 - 1 /p)-i/2 (1 + s) -Us 

Jo 

icciEo+atf^i+t)-^ 1 - 1 ^. 

Using (I4TT3D and (Oil) - ([4T25D . we obtain \v(-,t)\ H K {x) < C(E + ((t) 2 )(l+t)-^ . Combining 
this with (|4.25p . p = 2, rearranging, and recalling definition (|4.20|) . we obtain (|4.3|) . □ 

Proof of Theorem By short-time H existence theory, \\(y,ipt,ip x )\\ H K is continuous 
so long as it remains small, hence r\ remains continuous so long as it remains small. By 
(|4.3p . therefore, it follows by continuous induction that n(t) < 2Ct]q for t > 0, if rjo < l/4:C, 
yielding by (|4T20|) the result (fLT5]) for p = 2. Applying ()4T24"l) - (l435l) . we obtain (fTTT5l) for 
2 < p < for any < oo, with uniform constant C. Taking p* > 4 and estimating 

|Q| L2 , |i?| i2 , |5| L2 (f) < |(«,^,^)|i4<C^(l + t)-t 
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in place of the weaker (|4.22p . then applying Corollary 13.11 with q = 2, d = 1, we obtain 
finally 15j) for 2 < p < oo, by a computation similar (|4.24p ^ (|4.25p ; we omit the details of 
this final bootstrap argument. Estimate (I1.16P then follows using (13.3)) with q = d = 1, by 
(4.26) 

|^(*)|xp < CE + C((tf f (1 + t - + s yhds < C(l + t)^(JSb + C(t) 2 ), 

J o 

together with the fact that i) — u(x) = v(x — ip,t) + (u(x) — u(x — tp), so that \u(-,t) — u\ 
is controlled by the sum of \v\ and \u(x) — u(x — ip)\ ~ \if)\. This yields stability for 
[it — it|^i n ^K|t=o sufficiently small, as described in the final line of the theorem. □ 



5 Nonlinear stability in dimension two 

We now briefly sketch the extension to dimension d = 2. Given a solution u(x,t) of (jl.4p . 
define the nonlinear perturbation variable 

(5.1) v = u — u = u(x± + ij){x, t),X2,t) — u(xi), 
where 

(5.2) u(x,t) := u(xi + ip(x,t),t) 
and ip : M. d x K — > R is to be chosen later. 

Lemma 5.1. For u ; it as in (|5.2p . 

(5.3) u t + /'(u)^ - u ^xj = (dt - L) u'(xi)ii>(x, t) + ^ + 5 4 5 + T, 

3=1 i=i j=i 

where 

Rj = 0((\v,ifj t ,'ipx)\\(v,v x ,ipt,'4 , x)\), S := -vtp Xl = (\v\(\t/) x \), ^ := 0(\if) x \ 3 +\(v,ip x )\\ip XQ 

Proof. Similarly as in the proof of Lemma |4. II it follows by a straightforward computation 
Using the fact that ut + Ylj dfi(u)u Xj — . u XjXj = 0, it follows that 

Ut + d f j ( u ) u xj ~ Ux o x i = Uxiipt ~ u t tp Xl + df j (u)u Xl ^ Xj 

(5.4) ' ' _ j¥1 _ 

j¥i j 

where it is understood that derivatives of it appearing on the righthand side are evaluated 
at (x+tp(x,t),t). Moreover, by another direct calculation, using the fact that L(u'(x\)) = 
by translation invariance, we have 

(d t - L) u' (xi)i/j = u Xl ip t - u t i/j Xl +'^2df j (u)u Xl ifj Xj -y~)ux iXl il) X:j - ^(^i^jW 

m m 3 
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Subtracting, and using f|4.5j) and 
(5.5) 



"ay + V Xj = U X] + U xl lp Xj , 



u t + v t = u t + u xi ip t , 
so that 

~ — ( ( — | \ j 

^aij — V^^i ~r "sci / -. . i ; 

(5.6) ' 

u t -u t -v t = -(u Xl + UsJ , 

^ ~r Wx\ 

we obtain 

n * + ^2 df j (u)u Xj - ^2u XjXj = (d t - L)u'(xi)ip + v xi ip t - v t ip Xl 
i j 

+ ^2{df J (u)u xl - dp (u)u Xl )ip Xj 

~ ^(UxjX! ~ U XjXl )lp Xj - ^((Uxi - U Xl )ij Xj ) a 



Using v Xl ip t - v t ip Xl = (vipt)xi ~ (vip Xl ) t , 

dp(u)u Xl = f(u) Xl - dp (u)u Xi i/j Xi = f(u) Xl (l - ip. x ) - df '(u)u Xl ipl v 
and u XjX1 = (u Xj ) Xl - u X]Xl ip xl = (u Xj ) Xl (l - ip Xl ) + u X]Xx ^ 2 x ^ and rearranging, we obtain 

Ut + d f j ( u ) u xj ~ ^2 Ux i x i = ~ L )u'(xi)ip + (vipt)xi ~ (vip Xl ) t 
j j 

+ Y,U J (n)-P(u)) Xl )^ Xj 

~ ~ ^xi^Xj + ^2(u Xj ) Xl 1p Xl tp Xj 

+ ^u XjXl ip 2 Xl ^ Xj 

IpXilpx 



Noting that 

(P(u) - f j (u)) Xl )ip x . = ((P(u) - P(u)^ Xj ) Xl - {P(u) - P(u))i> XjXl , 
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fiu) x ^ xi ij) x . = (f(u)ip xi ip x .) xi - f(u)(ip Xl ip Xj ) xl , 

and 

with \f 3 (u) — f J {u)\ = 0{\v\) and \u Xj — u Xj \ = 0(\v\), we obtain the result □ 

Proof of Theorem The result of Lemma 15.11 is the only part of the analysis that differs 
essentially from that of the one-dimensional case. The cancellation and nonlinear damping 
arguments go through exactly as before to yield the analogs of Propositions 14.31 and (|4.2p . 
Likewise, we obtain a Duhamel representation analogous to (|4.18p ~ (|4.19p . forming a closed 
system in variables (v,t/j x ,t/jt)- 

To obtain the analog of Lemma 14.31 completing the proof of nonlinear stability, we can 
carry out a somewhat simpler argument than in the one-dimensional case, using Corollary 
13.11 with d = 2, q = 2 for all estimates, not only the final bootstrap argument, giving in 
place of (|4.24p the estimate 
(5.7) 

H;t)\ LP{x) < C(l + ty^-^Eo + CC(t) 2 [ (1 + 1 - sJ-CVa-Vp)^ _ s )~!(i + s )-^ ds 

Jo 

<c(£ + cw 2 xi + r (1 ~ 1/p) , 

\{^ t )(^t)\ Lv{x) <C{l+t)-^- l M-\E Q 
(5.8) + C((t) 2 [ (1 + t- 8 )-W*-V p \t - s)"5(l + s^ds 







<C( J E + C(i) 2 )(l + ^ (1_1/P) "5 
for divergence- form source terms, and 

v(;t)\ LP (x)<C((t) 2 I (l+t-sy^-^il + syUs 
Jo 

KCiEo + atf^l + ty^M, 



(5.9) H-M L ^<c«t?f Q 



K^,^)(-,t)u P(:C )<c(i+ir (1-1/p)- ^o 

(5.10) + C((t) 2 / (1 + t - s)- {1/2 - lM (t - s)-2(l + s )~2ds 

Jo 

< C(£ + CW 2 )(1 + 

for faster-decaying nondivergence-form source terms. 

We omit the details, which are entirely similar to, but substantially simpler than, those 
of the one-dimensional case. □ 
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